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ABSTRACT 

The dense molecular clump PI in the infrared dark cloud (IRDC) complex 
G28.34+0.06 harbors a massive protostellar cluster at its extreme youth. Our 
previous Submillimeter Array (SMA) observations revealed several jet-like CO 
outflows emanating from the protostars, indicative of intense accretion and po- 
tential interaction with ambient natal materials. Here we present the Expanded 
Very Large Array (EVLA) spectral line observations toward PI in the NH3 (J,K) 
= (1,1), (2,2), (3,3) lines, as well as H2O and class I CH3OH masers. Multiple 
NH3 transitions reveal the heated gas widely spread in the 1 pc clump. The 
temperature distribution is highly structured; the heated gas is offset from the 
protostars, and morphologically matches the outflows very well. Hot spots of 
spatially compact, spectrally broad NH3 (3,3) emission features are also found 
coincident with the outflows. A weak NH3 (3,3) maser is discovered at the in- 
terface between an outflow jet and the ambient gas. These flndings suggest that 
protostellar heating may not be effective in suppressing fragmentation during the 
formation of massive cores. 

Subject headings: ISM: individual (C28.34+0.06) — ISM: jets and outflows — 
ISM: magnetic flelds — stars: formation — stars: early-type — masers 

1. Introduction 

The majority of stars, from low to high mass ones are born in clusters (Lada & Lada 
2003). In a clustered environment, feedback from protostellar radiation and outflows can 
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play an important role in shaping their parental molecular cloud and thereby regulate fur- 
ther cloud fragmentation and star formation. In particular, protostellar feedback may be 
crucial for the formation of massive stars which typically contain masses at least one order 
of magnitude higher than the global Jeans mass in the molecular clump (Mj ^ 1 Mq for 
nn ^ 10^ cm"'^ and T ~ 15 K; Pillai et al. 2006; Rathborne et al. 2006). Calculations have 
shown that the accretion luminosity of a low mass protostar can heat gas to 100 K out to 
hundreds of AU, thereby raising the local thermal Jeans mass (Masunaga & Inutsuka 2000; 
Krumholz 2006). Subsequent numerical simulations indicate that protostellar heating from 
low mass stars in the cluster may change the equation of state (Larson 1985), and suppress 
the fragmentation (Krumholz et al. 2010), thereby aiding the formation of massive cores. On 
the other hand, protostellar outflows in the cluster can also inject momentum and energy 
into the local medium, and affect the evolution of the cloud (Quillen et al. 2005; Curtis et al. 
2010; Wang et al. 2010). 

To investigate the protostellar feedback, we observed the dense molecular clump PI 
in the IRDC complex G28.34+0.06. Located in the southern part of G28.34+0.06, PI is 
one of the dense molecular clumps among the filamentary network in the IRDC complex 
which stretches over 6 pc at a kinematic distance of 4.8 kpc (Carey et al. 2000; Rathborne 
et al. 2006; Wang et al. 2008). PI is a northeast-southwest (NE-SW) orientated molecular 
clump containing a large mass reservoir of 10^ Mq within a size of 1 pc. Wang et al. (2008) 
presented an NH3 mosaic of the entire IRDC, and found that the gas temperature decreases 
from 20 to 14 K toward the inner part of PI, indicative of external heating. A similar 
trend was found in the NH3 line width, suggesting a turbulence dissipation toward smaller 
scales. Our SMA observations resolved the PI clump into a string of five cores (named 
SMAl-5) containing 22 64 Mq in mass and ^0.1 pc in size; at higher angular resolution 
and better sensitivity, the cores are further resolved into ten compact condensations (named, 
e.g. SMA2a d) with masses of 1.4-10.6 Mq and sizes of ^0.01 pc, indicative of a hierarchical 
fragmentation (Zhang et al. 2009; Wang et al. 2011). In addition, all five cores are associated 
with highly coUimated CO (3-2) outflows. The large mass reservoir and the ongoing accretion 
indicate that PI will eventually form a massive star cluster after the accretion is completed. 
This young protocluster provides an opportunity to investigate how protostellar heating 
and outflows affect the parent molecular environment, and the formation of massive cores. 
The new EVLA observations presented in this Letter reveal the temperature proflles at 
unprecedented high angular resolution, allowing us to evaluate the protostellar feedback in 
detail. 
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2. 



Observations 



The NRAO^ EVLA was used to obtain the NH3 (J,K) = (1,1), (2,2), (3,3), as weU as 
H2O and CH3OH maser hne emission data in 2010 during the EVLA early science phase. 
Table 1 summarizes the observations. The data were calibrated in CASA^ and imaged using 
CAS A and Miriad^. 

The NH3 (1,1) and (2,2) observed in the EVLA C array has a channel width of 0.2 
km s~^. Wang et al. (2008) observed NH3 (1,1) and (2,2) toward the same phase center 
in the VLA D array with a channel spacing of 0.6 km s~^. We combined the calibrated 
visibilities of Wang et al. (2008) with the new C array data for imaging. To increase the 
signal-to-noise ratio in the final images, we averaged the channel width to 0.3 km s~^ and 
tapered the visibilities with a Gaussian function with a FWHM of 180 kA. We also made 
images with a 0.6 km s~^ channel width to compare with the images at 0.3 km s~^ resolution, 
and found that the latter does not introduce additional features. For the NH3 (3,3) images 
we also averaged the channel width to 0.3 km s~^ to improve the signal-to-noise ratio. 



Figure 1 presents moment 0, 1 and 2 maps of the NH3 (1,1), (2,2) and (3,3) main 
hyperfine emission. The NE-SW orientated filament is traced by the NH3 (1,1) and (2,2) 
emission extremely well. Ammonia cores coincide with all the SMA dust cores, indicating 
that NH3 traces the dense envelopes of the protostellar cores. As illustrated in the moment 
1 images of NH3 (1,1) and (2,2), there is a velocity gradient of 1.5 km s~^ over 1 pc along the 
filament. This gradient is, however, not uniform; it is enhanced toward southwest starting at 
SMA5. On the other hand, the velocity field from SMAl through SMA4 is quite uniform, at 
about 79 km s~-^. The moment 2 images show similar velocity dispersion^ of 0.4 0.8 km s~-^ 
along the filament as traced by NH3 (1,1) and (2,2). The velocity dispersion is enhanced to 
several local maxima of 0.6-0.8 km s~^ which coincide with all the SMA dust cores except 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

^ http://casa.nrao.edu 

^ http:/ /www. astro. umd.edu/~teuben/miriad, http:/ /www. cfa.harvard.edu/sma/miriad 

^For a Gaussian line profile, velocity dispersion a is related to FWHM line width AF as cr = AF/2/21n2. 



3. 



Results 



3.1. 



Ammonia Emission 
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SMA3. The enhancements may originate from a combination of spatially unresolved motions 
such as rotation, infall, and outflow. 

In addition to the main fllament, the NH3 (1,1) emission also reveals a fainter northwest- 
southeast (NW-SE) orientated fllament, crossing the main fllament at SMA5, making an X 
shape. This minor fllament is located at a slightly different velocity, around 82 km s~^. 
Compared to the main filament, the minor filament shows weaker emission in NH3 (1,1) and 
little emission in NH3 (2,2), indicative of even colder gas. The minor filament is evident 
in the previous NH3 (1,1) image at lower resolution (Wang et al. 2006; Zhang et al. 2009), 
however, it does not correspond to any dust emission feature at 350-1300 /xm (Rathborne 
et al. 2006, 2010; Zhang et al. 2009; Chen et al. 2010; Wang et al. 2011). Therefore, the 
minor filament is remarkably less dense than the main filament. Like the main filament, the 
minor filament consists of a velocity gradient of 2 km s~^ over 1 pc; it also hosts several NH3 
cores with sizes similar to the NH3 cores on the main filament. It has a velocity dispersion 
of < 0.4 km s~^ in general, but local enhancements of up to 0.5 km s~^ are also presented. 
We will discuss the nature of the minor filament in Section 4.3. In the following text, unless 
otherwise stated, we refer the main NE-SW filament as "the filament" . 

It is interesting to note that the clumpy NH3 (3,3) emission does not follow the dust 
filament at all. We visually identified five representative peaks from the (3,3) moment 
map, labeled as A E in Figure 1. None of these peaks coincide with any dust cores. In 
fact, the only one that lies on the filament is peak D which is between SMA4 and SMA5; 
other peaks are offset from the filament. Because NH3 (3,3) is excited at an energy level of 
125 K, these ammonia peaks trace local temperature enhancements relative to the generally 
cold molecular clump. Another important feature is the large velocity dispersion: all the 
peaks show centrally condensed distribution of cr > 0.8 km s~^, and up to 2.2 km s~^ in 
peak B. Peak B is located down stream of the red lobe of outfiow SMA2a, the most energetic 
outflow in PI (Wang et al. 2011). The centroid velocity on peak B is redshifted relative to 
the fllament, consistent with the redshifted CO emission. Indeed, all these NH3 (3,3) peaks 
seem to trace the footprint of the jet-hke CO (3-2) outflows revealed by the SMA (Wang 
et al. 2011). We will discuss this further in Section 4.1. 



3.2. Maser Emission 

We report maser detections based on the following criteria: (a) Peak flux >6(j; (b) 
Emission continues at least two channels with flux >4(j; and (c) Within the primary beam 
response of >0.3. We adopt these criteria because at some channels the images are con- 
taminated by the side lobes of a strong maser located at P2, well outside of the primary 
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beam. Using the criteria, we identify 4 water masers (Wl-4) and 1 methanol maser (Ml), 
as presented in Figure 2. Positions of the masers are obtained by fitting a two dimensional 
Gaussian to the brightness distribution of the maser emission. The positional uncertainty is 
<0'^02. All masers are either coincident with or close to infrared point sources at 8-24 fim. 

Comparing the maser positions and velocities with previous observations made with 
VLA (Wang et al. 2006, resolution -2'') and GBT (Chambers et al. 2009, resolution -^30''), 
we find that Wl, W3, and W4 are new features, while W2 and Ml are spatially and spectrally 
coincident with two masers reported before. Both W2 and Ml have been varying in peak 
fiux during the last several years: W2 was 2.4 Jy on 2005 September 23/24, then brightened 
to 6 Jy on 2006 December 5, and dimmed dramatically to 0.03 Jy on 2010 November 24; Ml 
dimmed from 0.24 Jy on 2006 December 5 to 0.02 Jy on 2010 November 24. On the other 
hand, Wl and W3 were either newly emerged or experienced brightening between 2006 and 
2010, otherwise they should have been detected by Chambers et al. (2009). W4 is below the 
detection hmit of Chambers et al. (2009). 

Among all the five detected masers, only two (Wl, W2) are associated with NH3 emission 
above 6a significance. Wl and W2 locate right on top of the dust condensations SMA2a and 
SMA2b, respectively, and thus are presumably excited by the protostars embedded in the 
same dust condensations. The position discrepancy between the H2O maser with its exciting 
protostar is within a fraction of 1000 AU. 

4. Discussion 

4.1. Outflow Heating 

Intensity ratios of the main and hyperfine components from multiple NH3 (J,K) transi- 
tions constrain rotational temperature which well approximates kinetic temperature in the 
regime of <20 K (Ho & Townes 1983; Walmsley & Ungerechts 1983). In the following analy- 
sis, we assume that the rotational temperature is the same as the kinetic temperature when 
the former is <20 K. Following the procedure of Ho & Townes (1983), we compute a rota- 
tional temperature (T21) map based on the NH3 (1,1) and (2,2) images. The NH3 images 
were first smoothed to a common beam of 2''. 8 before computing the T21 map. The uncer- 
tainty of T21 is estimated to be 3 K. Figure 3(a) shows T21 with the dust condensations and 
the CO (3 2) outflows superposed. The T21 map reveals a highly structured temperature 
profile ranging from 8 to 30 K, opposite to the smooth temperature profiles observed in other 
IRDCs (Ragan et al. 2011). This may be due to beam dilution since the data in Ragan et al. 
(2011) have a larger beam. A number of remarkable enhancements {T21 > 26 K) are present 
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among the entire map, but surprisingly none of the temperature enhancements coincide with 
any dust condensations. Instead, all the T21 enhancements spatially match the footprint of 
the outflows. 

Since the NH3 (3,3) is more widely spread than the (1,1) toward the northwest where 
the signal-to-noise ratio of (1,1) prohibits us from computing opacity, it is difficult to derive a 
reliable rotational temperature map from the NH3 (1,1) and (3,3) lines. Instead, we compute 
the intensity ratio of the main hyperfine of (3,3)/(l,l), i?3i, as plotted in Figure 3(b). To 
compute the intensity ratio, the NH3 (1,1) image was first smoothed to the beam of the 
(3,3) image. The intensity ratio is positively correlated with the rotational temperature (Ho 
& Townes 1983; Zhang et al. 2002). For reference, assuming optically thin emission and 
the same abundance for the ortho and para NH3 (Zhang et al. 2007), Rsi = 0.5, 2, and 5 
correspond to rotational temperatures of T31 = 31, 55, and 109 K, respectively. Like the 
T21 map, the R^i map is also highly structured. R^i varies continuously from 0.3 to 12.5 
over the map. Again several hot spots represented by Rsi peaks (corresponding to NH3 (3,3) 
emission peaks A-E) match the outfiow footprint and are offset from the dust condensations. 
In addition to the peaks, there is a low Rsi valley along SMA2 to SMA3. Optically thin 
thermal NH3 (3,3) emission should have Rsi ^ 6.4 (Ho & Townes 1983). Higher R^i values 
are seen toward peaks A-C, indicative of non-thermal emission. Toward peak B, the (3,3) 
emission is weakly amplified (Figure 4), resembhng an NH3 (3,3) maser similar to those 
found in other massive star formation regions where outflows are present (e.g. Mangum & 
Wootten 1994; Zhang & Ho 1995; Beuther et al. 2007; Brogan et al. 2011). However the 
spatial resolution is insufficient to constrain the brightness of this emission feature. 

The morphological match between the heated gas and the outfiow footprint suggests 
that the heating in the PI region is due to protostellar outfiows. Indeed, all the hot spots 
exclusively coincide with large NH3 (3,3) line widths compared to relatively small line width 
toward the dust condensations (Figure 4), likely indicative of gas motion and extra turbulence 
injected by the outfiow shocks. The outfiow heating has two important features. First, the 
heating alters thermal property of the gas throughout the clump, from a distance 0.1 pc up 
to 0.5 pc away from the driving protostars. Second, the heating is localized and far from 
well mixed. The temperature at the low Rsi valley is not high enough to suppress further 
fragmentation into Jeans mass of 1 Mq (Wang et al. 2011). We probably have witnessed 
through EVLA and SMA the very beginning of an outfiow induced gas heating. 



4.2. Radiative Heating 



It is apparent that protostellar outflows can introduce signiflcant heating to the ambient 
gas. What is the effect of protosteUar radiation? Is radiative heating from low mass proto- 
stars suflacient to suppress fragmentation and aid the formation of massive cores? The fact 
that all dense cores correspond with low Rsi is intriguing. In addition, the temperature map 
in Figure 3(a) indicates that dense cores have typical temperatures of 18 K, and are offset 
from the high temperature peaks. The temperature distributions suggest that the heating of 
low mass protostars may not play a significant role in the formation of dense cores detected 
with the SMA. Otherwise, massive cores should be found at the peaks of the temperature 
distribution. Longmore et al. (2011) made a similar statement for a hot molecular core 
G8.68-0.37. 

It is also interesting to note that the protostars embedded in the SMA cores have not 
produced significant heating. SMAl to SMA4 have a rather uniform T21 of 18 K, while SMA5 
has a lower T21 of 9 K. We caution that the T21 map has a 2'^8 beam, so any temperature 
structure has been smoothed to a spatial scale of 0.065 pc. How much luminosity is required 
to heat the gas to 18 or 9 K at a radius of 0.065 pc from the central source? Assuming 
thermal equilibrium between gas and dust, the luminosity can be estimated following the 
equation (Scoville & Kwan 1976; Garay & Lizano 1999): 



where /3 = 1.5 is the dust emissivity at far infrared wavelengths (Rathborne et al. 2010), 
/ = 0.08 cm^g"-^ is the dust emissivity at 50 /xm (Scoville & Kwan 1976), and r = 0.065 
pc is the radius from the central protostar. We estimate a luminosity of 30 Lq for each 
protostar embedded in SMAl, SMA2a, SMA3, and SMA4a, and 1 Lq for the protostar 
SMA5. These values should be taken with caution because the beam dilution may introduce 
large uncertainty in both temperature and radius. Nevertheless, we note that the luminosity 
of SMA2a is in rough agreement with the luminosity inferred from SED for the 24 fim source 
in the SMA2 vicinity (^10^ L©, Wang et al. 2008). A lack of significant heating from the 
growing protostars is consistent with non detection of line emission from organic molecules 
(Zhang et al. 2009; Wang et al. 2011). 



The origin of the low-density minor filament is rather interesting. A low-density streak 
perpendicular to a high-density oblate region is expected if magnetic field is dynamically 




4.3. Origin of the Minor Filament 
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important compared to self-gravity and turbulence during the formation of the system (Li 
et al. 2011). In this case, gravitational contraction along the field lines is more efl&cient, 
resulting in a high-density region perpendicular to the magnetic field, while on the other 
hand turbulence is also channeled by the field and gives rise to streaks aligned with the 
magnetic field in the lower-density region. This phenomenon has been seen in simulations 
(Price & Bate 2008; Nakamura & Li 2008) and in observations of NGC 2024 in Orion (Li 
et al. 2011). In G28.34-P1, the X-shaped filamentary system projected on the sky is almost 
perpendicular. If the system were formed from the same mechanism, one should expect a 
mean magnetic field along the minor filament and perpendicular to the main filament. Based 
on the JCMT dust polarization data cataloged by Matthews et al. (2009), we derived a mean 
B-field direction of P. A. = —25°. 6 ± 45° as labeled on Figure 1. Because the main filament 
has a P. A. of 48° (Wang et al. 2011), the separation angle between the mean B-field and 
the main filament is 73°. 6. The mean magnetic field at the clump scale is indeed along with 
the minor filament and roughly perpendicular to the main filament. The large deviation 
in the B-field direction may indicate that the strong outfiows have started to disturb the 
initially ordered field fines in the inner part of PI (Li et al. 2009). Therefore we speculate 
that the minor filament originates from turbulence anisotropy during the formation of the 
X-shaped filamentary system, and that the formation process has been governed by the 
interplay between strong magnetic field, self-gravity, and turbulence. 

We thank Edward T. Chambers for providing the spectra of the masers detected in 
G28.34+0.06. K.W. acknowledges the support from the SMA predoctoral fellowship and the 
China Scholarship Council. Q.Z. acknowledges the support from the Smithsonian Institution 
Endowment Funds. This research was supported in part by the NSFC grants 10733030, 
10873019 and 11073003. 
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Table 1. Observational Parameters 



Parameter /Line 




NH3 (1,1) 


NH3 (2,2) 


NH3 (3,3) 


H2O/CH3OH 


Observing date 




2010Oct09 


2010Octll 


2010May09 


2010Nov24 


Configuration 




C 


C 


D 


C 


Rest frequency (GHz) 




23.694 


23.722 


23.870 


22.235/24.959 


Primary beam {') 




2 


2 


2 


2 


Gain calibrator 




J1851+005 


J1851+005 


J1851+005 


J1851+005 


Flux calibrator ^ 




3C48 


3C48 


3C48 


3C48 


Bandpass calibrator 




3C454.3 


3C454.3 


3C454.3 


3C454.3 


Phase center ^ 




I 


I 


II 


II 


Integration time (minute) 




95 


95 


21 


15/15 


Bandwidth (MHz) 




4 


4 


4 


4 


Polarization 




dual 


dual 


dual 


fuU 


Naive channel width (km s" 




0.2 


0.2 


0.2 


0.8 


Final channel width (km s~ 




0.3 


0.3 


0.3 


0.8 


Synthesized beam ('') ^ 




2.7 X 2.2 


2.8 X 2.0 


4.8 X 2.7 


1.3 X 1.0/1.5 X 0.9 


RMS noise (mJy/beam) ^ 




1.2-1.5 


1.2-1.5 


3 


1.2-1.5/0.8 


Convention factor (K/Jy) 




367 


388 


165 


1939/1484 


Weighting method 




natural 


natural 


natural 


robust / natural 



^Absolute flux is accurate to about 5%. 

^Phase center I at (R.A., decl.)j2ooo = 18^42^50^ -04°03'30''; H at 
18M2"^50^82, -04°03al'^3. 

^Final NH3 (1,1) and (2,2) images are made by adding archival D array data (Section 2). 

^la rms noise measured in final images. 
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NHg( 1,1) NH3(2,2) NH3(3,3) 




18''42'52" 5l" 50" 49" 
a (2000) 



Fig. 1. — Moment maps of the main hyperfine emission of NH3 (1,1), (2,2) and (3,3), shown 
on the left, middle and right columns, respectively. The contours are plotted in steps of 6(7, 
which is 0.5 mJy beam"-*^ km s"-*^ for (1,1) and (2,2), while 1 mJy beam"-*^ km s~-^ for (3,3). 
Units of the color bars are km s~^. The stars mark the five dust cores (SMA1^5) which are 
further resolved into smaller condensations (the black contours in Figure 3). The blue line 
on the bottom-left panel indicates the direction of the mean magnetic field (Section 4.3). 
The filled ellipses on the bottom-left corners represent synthesized beams. 
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Fig. 2. — Detected water (Wl-4) and methanol (Ml) masers. Left: Spitzer composite image 
(red/green/blue = 70/24/8 /xm) showing the IRDC G28.34+0.06 complex overlaid with the 
locations of the masers. The Spitzer data are taken from the GLIMPSE and MIPSGAL 
legacy projects (Churchwell et al. 2009; Carey et al. 2009). Right: Spectra of the masers. 
Brightness temperatures are corrected for primary beam response. 
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J2000 Right Ascension 




J2000 Right Ascension 



Fig. 3. — Morphological match between heated gas and protostellar outflows. Plotted in 
color scales are (a) T21, the rotational temperature derived from NH3 (1,1) and (2,2), and 
(b) i?3i, the integrated intensity ratio of NH3 (3,3)/(l,l). The arrows and blue/red contours 
outline the blue-/redshifted CO (3-2) outflows, and the compact black contours represent 
the flve groups of dust condensations, SMAl-5, from upper-left to bottom-right (Wang et al. 
2011). The two white crosses on top of SMA2a and SMA2b represent water masers Wl and 
W2, respectively. On the bottom-left corners are the beams for dust continuum, CO (3^2), 
and T21 (i?3i), from the smallest to the largest respectively. Note that masks of 5a were 
applied to both maps, and additional manual clipping was also performed to Rsi in order to 
avoid noisy edge pixels. 
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Fig. 4.^ EVLA NH3 spectra extracted from SMA dust peaks (SMAl-5) and NH3 (3,3) 
peaks (AD). Note at peak B a spectroscopically unresolved NH3 (3,3) maser at about 84 
km s~^. (Panel E was truncated in formal publication due to the space limit of ApJ Letters.) 
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